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Abstract. For a large class of linear neutral type systems the problem of 
eigenvalues and eigenvectors assignment is investigated, i.e. finding the system 
which has the given spectrum and almost all, in some sense, eigenvectors. 



1. Introduction 



One of central problems in control theory is the spectral assignment problem. 
This question is well investigated for linear finite dimensional systems. It is impor- 
tant to emphasize that the assignment of eigenvalues is not sufficient in several cases. 
,—h ■ One needs also the assignment of eigenvectors or of the geometric eigenstructure. 

For infinite dimensional problems (delay systems, partial derivative equations) the 
problem is much more complicated. 

Our purpose is to investigate this kind of problems for a large class of neutral 
type systems given by the equation 

(f.I) z(t) = A-!z(t-l)+ A 2 (e)z(t + 6)d0+ A 3 {6)z(t- 



where z(t) G M. n and A-\, A2, A3 are n X n matrices. The elements of A 2 and A3 



_±. ' taking values in L 2 {— 1, 0). The neutral type term A—\z(t — 1) consists on a simple 

delay, while the other include as multiple as distributed delays. The behavior of 
such systems can be described mainly by the algebraic and geometric properties of 

^f\ . the spectrum of the matrix A_\ (cf. [UED- 

It is well known that then spectral properties of this system are described by the 
characteristic matrix A(A) given by 

\e X0 A 2 {0)d9- / e xe A 3 (9)d9. 



S : 

The eigenvalues are roots of the equation det A (A) = 0. The eigenvectors of the 
system (more precisely of the functional operator model of the system) are expressed 
through the matrix A(Afc), where Afc is an eigenvalue. In fact the problem of an 
assignment of an infinite number of eigenvalues and eigenvectors is reduced to a 
problem of assignment of singular values and degenerating vectors of an entire 
matrix value function A(A). It is remarkable [5] that the roots of det A(A) = 
are quadratically close to a fixed set of complex number which are the logarithm 
of eigenvalues of the matrix A_i. Moreover, the degenerating vectors of A(Afc) arc 
also quadratically close to the eigenvectors of the matrix A-\. 
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In this paper we investigate an inverse problem: 

What conditions must satisfy a sequence of complex numbers {A} and a sequence 
of vectors {v} in order to be a sequence of roots of the characteristic equation 
det A(A) = and a sequence of degenerating vectors of the characteristic matrix 
A(A) of equation (|1.1[) respectively for some choice of matrices A_j, A 2 {6), A 3 (9) ? 

One of the possible application of this problem is to investigate a vector moment 
problem via the solution of the exact controllability property for a corresponding 
neutral type system by a relation devlopped in |3]. 

The present paper is a detailed version of the short note published in Comptcs 
Rcndus Mathcmatiques [7]. 

2. Operator form of perturbation 
We consider neutral type systems of the form 

(2.1) z(t) = A- t z{t - 1) + / A 2 (6)z(t + 9)d6 + A 3 (6)z(t + 



where A-\ is a constant n x n matrix, A 2 ,A 3 are nx n matrices whose elements 
belong to L 2 [-1,0]. 

As it is shown in [3], [5] this system can be rewritten in the operator form 

di{z t (-)) =A {z t (-)J' 

where A : D(A) -> M 2 = C" x L 2 ([-1,0],C"), 

D(A) = | ( * ) | tp(-) e ^([-l.Oj.C 1 ), y = <p{0) - A_i^(-1)} c M 2l 
and the operator A is given by formula 

( y \ fJ Q _ 1 A 2 (9)ip(e)de + j Q 1 A 3 (9Me)de\ 

W-)J I &(•) )■ 

This operator is noted A instead of A if A 2 (9) = A 3 (0) = 0. The operator A is 
defined on the same domain D(A). One can consider that the state operator A is 
a perturbation of the operator A, namely 

y \ =1 ( y \^ fj^A 2 (6Me)d6 + j° 1 A 3 (6M6)de 



M-)J \<p( 

Let B : C" -> M 2 be given by 

'2/ 



and Vo :D(A)^ C™ by 

(2.2) Vo {j(-)) = I A ^ e ^ d9 + f MdM0)d6. 

Then .4 = .4 + BqVq. Denote by Xj± the set D(A) endowed with the graph norm. 
Let us show that Vq browses the set of all linear bounded operators C(Xj\, C n ) as 
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A 2 (-), As(-) run over the set of n x n matrices with components from L/2[— 1,0]. 
Indeed, an arbitrary linear operator Q from £(Xa,C") can be presented as 

o( (.)) = Qiv + QM-) 

= Q 1 (<p(0)-A- 1 <p(-l)) + J A 2 {p)<p{0)&0 + J A 3 (d)cp(6)d6, 

where A 2 {-), A 3 (-) are (n x n)-matrices with component from L 2 [—\, 0] and Qi is a 
(n x n) matrix. Let us observe that 

/0 /-O 

^j(#)d0 + / ¥>(0)d0, 

/-O 



<p(0) = / (e + l)^(9)d6+ / y>(0)d0 



and denote 

A 2 {6) = A 2 (6) + (6+l)Q 1 -eQ 1 A^ 1 , 
A 3 (6) = A z {e) + Q 1 -Q 1 A_ l . 

Then, with these notations, the operator Q may be written as 

0( o)=/ A 2 (8)^(9)d8 + J A 3 (0)<p(0)d6. 
Hence formula (|2.2|) describes all the operators from £(X^,C n ). 

3. AN EQUATION FOR EIGENVALUES AND EIGENVECTORS OF THE 
CHARACTERISTIC MATRIX (SPECTRAL EQUATION) 

Consider the operator A = A + BqPq and assume that Ao is an eigenvalue of A 
and xq is a corresponding eigenvector, i.e. 

(3.1) (A + BqVo)x = Xox . 

Let us assume further that Ao docs not belong to spectrum of A and denote by 
R(A, Ao) = (A — Ao/) -1 , with this notation (|3.ip reads as 

(3.2) xo + R{A,\o)B o Vox a = 0. 

Let us notice that vq = VqXq ^ 0, because Ao ^ a {A). Then applying operator Vo 
to the left hand side of (13.21) wc get 



v + VoR(A,Xo)B o v o = 

This equality means that Ao is a point of singularity of the matrix- valued function 
F(X) = I + VqR(A, A)£>o and vo is a vector degenerating F(Ao) from the right. 
Let Wq be a non-zero row such that 

(3.3) w* F(X ) = 0. 

In order to describe the vector wq, let us find first another form for the matrix 
F(X). For any v G C™ we denote 

R(A,X)B v= ( y V X{o(A), 



Then 



This gives 
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\«) -*•-(!! 



„ = -A,, ^|> - W = 0, 



and, as ( ¥%) £ D(A), we obtain j/ = tp(0) — A-i<p(— 1). Therefore 

99(6*) = L>c Ae , £> e C", 

and 

v = -A(/-c" A J 4_i)L». 

Since the matrix (/ — c~ A j4_i) is invertible (A ^ a(A)), then we get 

<p(e)=-^-(I-e- x A_ 1 )- 1 v, 
and hence 

/ 

This formula and (|2.2j) implies 
w + 7> i?(A A)S w = 

= v-(f c xe A 2 (6)d8- J ^-A 3 (6)de) (I - e^A^v 

and hence 

where 

/o /.O 

Ac Ae A 2 (6»)d6i- / c Ae A 3 (0)d0 

is the characteristic matrix of equation (jl.ll) . Thus, the equality (J3.3I) for Ao £■ c(.A) 
is equivalent to 

Wo *A(Ao)=0. 
Summarizing we obtain the following 

Proposition 3.1. Let Ao do not belong to a (A). Then the pair (Ao,u»o), u>o £ 
C n , wo ^0, satisfies equation (|3.3[) i/ and on??/ i/ Ao is a root of the characteristic 
equation 

detA(A) =0 
and wo* is a row-vector degenerating A(Ao) from the left, i.e. wo*A(\q) = 0. 

Thus, one can consider the equation w*F(X) = as an equation whose roots 
(A, w) = (Ao, wo) describe all eigenvalues and (right) eigenvectors of the character- 
istic matrix A(A). 
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4. A COMPONENT-WISE REPRESENTATION OF SPECTRAL EQUATION 

We recall spectral properties of operators A and A* obtained in [U \5\ . We will 
assume that the matrix A^i has a simple non-zero eigenvalues \x\,Hi, /13, . . . ,fi n . 
In this case the spectrum a(A) consists of simple eigenvalues which we denote by 

A™ = ln|^ m | + i(Arg / u m + 2nk), m = 1,. . . ,n, k E Z, 

and of eigenvalue Ao = 0. First assume [x m ^ 1, m = 1, ...,n. Then the corre- 
sponding to A™ eigenvectors are of the form 

(4.1) vZ=(~xZ V fceZ, m = l,...,n, 

where y\ , . . . , y n are eigenvectors of A_\ corresponding to /j,i, fJ,2, ■ ■ ■ , H>n- The eigen- 
spacc corresponding to Ao = is n dimensional and its basis is 

(4.2) ^( (1 ^' )yj )' i = l>-,»- 

If some /j, m say /ii equals 1, then Ap = Ao = 0. In that case the eigenspace cor- 
responding to is (n + l)-dimcnsional and its basis consists of n eigenvectors 
ip®, j = 1, . . . , n, given by (|4.2|) . and one rootvector 

All the vectors iffi given by l|4.1|) . except of ^5, are still eigenvectors of .4 and 
correspond to eigenvalues A™. In both cases the family $ = {^™} U {£>?} forms a 

Riesz basis in the space Mi. Denote by \& = < ^™ > U < V? f the bi-orthogonal basis 
to $. Then 

(z /A" 
"L k j , m = l,...,n, fceZ, A ^0, 
e- A ™ e z n 

where A is the complex conjugate of A and 

* - (;■ 

if Ag = 0, here z m are eigenvectors of matrix A*_ x such that 

It is easy to see that if)™ are eigenvectors of the operator A* corresponding to 

eigenvalues A™. 

Our nearest goal is to rewrite the matrix F(Xq) in the basis $. Let us first 
rewrite the expression R(A, Ao)£>o- We present Bq as a matrix Bo = (&?, b®, ■ ■ ■ , b„) 
with infinite columns b® which are vectors from Mi of the form 

b\ = f jj, i = l,...,n, 
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where ei is the canonical basis of C™. Then 

&? = £(^;r) fir +£(*?,#)$, 

k,m j 

and in the case \x m ^ 1, m = 1, . . . , n, we have 



i = l,...,n, 



k,rn k j 



'^, t = l,...,n. 



Taking into account the form of eigenvectors f|4.3[) we find that for all to = 1, . . . , \ 
bit 



k , 



6^^r)) = ~^((ei,^),...,(e„, Zm )) = ^/Ar. 
A,. 



So we obtain 

(4.4) 



^^^E^^^ 4 ^^^'^''--'^^) 



fc,m A k ~ A A fc 



if all the numbers A™ ^ 0. If /ii = 1 and Aj = then in the sum (|4.4[) the term 
corresponding to k = 1, to = 1, is replaced by 



1 



-<Ad 



1 



\A — Ao (A — Ao) 2 

In what follows, we shall use the notation 



~o 1 * 



ftk 



A?, Aj?^o, 

i, xt = o, 



m = 1, . . . , n, fcgZ. 



Let us observe that the expression w^Vo is a linear bounded functional on the space 
Xj±, i.e. WqVq G £(Xa, C). The representation of u>g in the basis Zj is as follows: 



< = ^2 a 3 z h 



Consider now n functionals z*Vq G £(Xa,C), j = 1, ...,n. One can decompose 
them in the basis ^>: 



k.m 



where 

(4.5) 



E 



pL 



& 



< 00, 



771 = 1, . 



In the sequel, we shall assume p^ = 0, i = 1, . . . , n. This means that from now we 
consider perturbations Vo satisfying the condition 



(4.6) 

Then we have 



II: 



A 3 {0)d6 = 0. 



*v = e «i*;p<> - E «>EW 
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From this relation and expression (|4.4j) we obtain the equality 

m \ k Pk A fc ~ A j 



fco,m 



a,- > — — x - z 

j k ,m Pk A fc ~ A 

With these notations, the equation (|3.3[) reads 



Pfe.m 1 
a,— — x — z„ 



= w* (I + V n(A, Xo)Bo) = E «m< + E E ^ , 

m m \ i,fc Pfe A fc - A 

Thus, the condition for a pair (Ao,u>o) to satisfy the spectral equation can be 
rewritten in the form of the following system of n equations: 

(4.7) a m = - y^ a 3 \ -^- x , m = 1, ...,n, 

fcez V Pk A k ~ A 0/ 

j=l,...,n 

where for any fixed couple m,j, the needed n-tuple Ipj, m > satisfies (14.51) . 

5. Conditions for spectral assignment 

Now we discuss the following question: 

What conditions must satisfy a sequence of complex numbers {A} and a sequence 
of vectors {v} in order to be a sequence of roots of the characteristic equation 
dot A(A) = and a sequence of degenerating vectors of the characteristic matrix 
A(A) of equation (|1.1[) respectively for some choice of matrices A_i, A 2 (9), A 3 (8)? 

We will assume that the corresponding operator A has simple eigenvalues only. 
Let us remember that we assumed earlier that all eigenvalues of matrix A-\ arc 
also simple. Then one can enumerate those eigenvalues as {A™} U {A°} , m,j — 
1, . . . ,n; k £ Z, where (see [5J Theorem 1]) the sequence {A™} satisfies 



(5.i) e| a ™- a ;' 

fc,m 



< 00. 



Denote by {y™} U {v?} > m >3 ' = lj • • • 7 1V , k £ Z, corresponding eigenvectors of A. 
Then (see [5J Lemma 13, Theorem 15]) these vectors form a Ricsz basis in M 2 which 
is quadratically close to the basis {<p™} U {£>?} if we assume that the corresponding 
elements have the same norm. Eigenvectors {(/?} has the form 



K A_ lV 

ip= ■ 



with A(A)u = 0. Therefore the fact that the basis {ip} and {(/?} are quadratically 
close implies the condition 

(5.2) YjIK 1 -Vnif <oo, m = l,...,n, 
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where A(A™)v™ = 0, m = 1, . . . , n; k G Z, and y m are eigenvectors of A-\ corre- 
sponding to fi m as in (|4.1j) . If we apply the same arguments for the dual system 

/o rO 

A* 2 (6)z(t + 8)d6+ / A* 3 (6)z(t + 6)d6, 

then we obtain the symmetric condition 

(5.4) ^||< l -z m || 2 <oo, m=l,...,n, 

k 

where (w™)*A(A^) = 0, m = 1, . . . , n; fc G Z. 

Our further goal is to show that conditions (|5.1I) , (|5.4[) or (j5.ll) , (|5.2[) arc almost 
sufficient for couples of sequences {A} , {w} or {A} , {v} to be spectral ones for the 
system (jl.ll) . 

Let us consider a sequence of different complex numbers {A™ } , mo = 1, . . . , n; 
&o G Z, satisfying (|5.1j) . We also assume that the index numbering of {A} is such 
that if A™ = AJ™ then m = m, fc = fc- To begin with, however, wc put AJ™° ^ A^ n , 
for all fc, fco G Z, ?7i, ??io = 1, . . . , n. 

Let now {A™ } be simple eigenvalues of operator A = A + BqVo, where Vo is 
given by (f2^|) . in which matrix A 3 (8) satisfies (|4T6|) . Then A(A^°) = and let 
{ui™° } be a sequence of the left degenerating vectors of A (A™ ), i.e. 

K7)*A(A-') = 0. 

We assume that the sequence w™° satisfies (|5.4|) . For all indices too = 1, . . . , n; 
fco G Z, consider decompositions 

n 

(«cr = E°iM- 

Then condition ()5.4|) is equivalent to 

(5.5) ^ \ a mm | < 00) m ¥" "to, ^ \ a n"rn ~ A < °°> TO, TO = 1, . . . , 77. 

k k 

Let us rewrite now relations (|4.7[) for A = A™ and u>o = w™° . 

We now consider the space £2 of infinite sequences (columns) indexed as {afc} feGZ 
with a scalar product defined by ({a^} , {b^}} = Ylk a kbk- From the relation (|4.5j) 
we obtain that vectors 

( — j ^ 

•j I I k.m I 

Fm = - 



— m 

/3 fc 



fcEZ, j, m— l,...,n 

belong to £2- One can also easily see that 
f 1 1 



A fc A k a ) feez 



G ^2, m, mo = 1, • • • , n; ko G Z. 



Then, putting Ao = A™ and wq = w™° in the equations (|4.7[) . we obtain 

(5-6) a^ jmo = ^ a^ ( <^ ~^ — - \ ,p° m ) , ™, to = 1, . . . , n; k G Z. 



\ 7TT, \ m 
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Now we would like to rewrite relation ()5.6|) in a vector-matrix abstract form. In 
order to do that, we introduce a more convenient notation. Denote 



n, 






m, m = 1, . . . , n. By S* mmo = {s"™° } fe feoGZ , m, m 
nite matrices with elements 



,n. 



we denote infi- 



mm 
b k k 



1 



*™ - K Q ° 



and by A,- mo , j, mo — 1, . . . , 71, inhnitc diagonal matrices 



^jmo — < ^ 1 ' 1 



ia s{ a X} 



fcoez 



With these notations relations (15.61) can be rewritten as 

(5.7) 



a r , 



/ j Ajrna £>mm Pri 
J = l 

m, mo = 1, . . . , n. Now let us fix index m and consider n equations (|5.T[) with this 
index and mo = 1,2, . . . , n. Consider another infinite diagonal matrix 



A m = diag{Ar-Ar} 



A-e: 



and multiply both sides of the m-th equality (|5.7I) (for mo = m) by this matrix 
from the left. This gives the following system of equalities 



(5.8) 



/ j Ajmo bmmoPn 



mo = 1, . . . , n, tuq ^ m, 



i=i 






^m^mm: 



where we used the fact that diagonal matrices commute : A m Aj m = Aj m A m . 
Finally, we introduce block matrix operators 

j4ll<Sml . . . A m i S ml ■ ■ ■ A n iS m l 



Dm = 



^-lm^-ni^n 



A A ^ 



mm 1 ^m L - f mm 



-ft-ln^mn • • • ^-mn^mn 

and present (|5.8|) in the form 



D„ 



A A <? 



^-nn^ran 



m = 1, . . . , n, 



«ml 



Let us observe that both vectors (p m , 
long to £rj = ^2 x ^2 x ... x £ 2 (see 






p^) r and (a m i . . . A m a mm . . . a mn ) be- 
. (|5~5j0 . Therefore the system (j5~5| is 
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solvable if and only if the vector (a m i . . . A m a mm . . . a m „) T belongs to the image 
of operator D m as an operator from £% to £%■ ^ n the sequel we show that, for all 
sequences {A} satisfying. (|5.ip and for almost all sequences {w} satisfying (|5.4p . 
operators D m , to = 1, . . . , n, are bounded and with bounded inverse operators from 
£2 to £%■ This means that the spectral assignment problem is solvable. In the 
further argument we use the following 

Proposition 5.1. Let {Afc} fceZ be a sequence such that 

y^ \Xk - a + i(& + 27rfc)| 2 < 00, 
fcez 

for some a, b £ R. Then the family x^ ht ) kl - 7 forms a Riesz basis in £2(0, 1). 

There are several ways to prove this classic result (see [1]). It may be be obtained, 
for example, from the Paley- Wiener theorem [2] and Lemma II. 4. 11 pQ. 
Next we prove the following preliminary result. 

Lemma 5.2. 1. For to 7^ m operators S mmo are bounded as operators from £(£2) 
and have bounded inverses. 2. A m S mm is a bounded operator from £(£2) and has 
a bounded inverse. 

Proof. Let {(fk} , {(f>k] , k £ Z, be two Riesz basis of a Hilbcrt space H and let R 
be a bounded operator with a bounded inverse, such that R(fk — <fiki k £ Z. 
For / £ H we have 



Then 



/ = X) a i W ' R f = ^2 a i Ri Pi 



Ripj =ipj = X c jk ip k = X (tPj^k) <Pk, 3 £ Z, 



where {V'fclfcez ^ s ^ ne bi-orthogonal with respect basis to {^fclfccz- Hence 
Rf = X a i X &i ' ^ ^ k = X) X a i &i > V>fc) ¥?fc = X bklfk ' 



where b k = £\ a 3 (<Pj,ipk) ■ 

This means that the infinite matrix R corresponding to R in the basis {</>&} is 
of the form 



where 



R 



R= {hj = (fj^k)} kez , 

J6Z 



, { aj }, {bj} £ £ 2 - 









a_i 




b-i 


a 


= 


bo 


a\ 




h 









Let now H = 1*2(0, 1) and {^fc} feGZ be a Riesz basis of the form tpk = c k t for 
some to = 1, ...,n. Let now {fk} be a Riesz basis which is bi-orthogonal to 
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< ipk = e A fc * > for some too = 1, . . . , n (the fact that {ipk}k<£Z i s a Riesz basis 
of £2(0, 1) follows from Proposition ^. ip . One has 



r° ~ 1 / 



J AC 



e A , " A * - 1 



Ar-Ar 



^m6 



,-ai 



i.e. -R = {rfcj} fceZj where 

J'GZ 



%,■ = Sfej Ut m e A <= - l) , k, j E Z. 



Thus 



-ft- ^mrriQ ^mmo 3 



where s m7nQ is the infinite matrix 









: Mme - 1 - 


- 1 





: 




/J m e" A » 


: 

















/i m e Al - 1 



1 > is bounded and separated 



Hence, we have the following alternative: 

1. If mo 7^ m, then the sequence < u m e _A i- 

L J fcGZ, 

from 0, i.e. e mmo : ^2 — ► ^2 is a bounded opeartor with a bounded inverse. Hence, 

(."■") "rotnj = e mm "- 

2. If to = toq, then n m & Xk —>!,&—> 00, moreover 



MmC A ™ - 1 



e A ™~ A ™ - 1 



f ur-Ar 



AJT - A£ 



= Af-Af l + o Af-A 



Therefore, 

where Q m = diag ( 1 + o(A™ — A™) J has a bounded inverse, so 

(5.10) A m S mm = Q-^ft 

From (|5.9[) . (|5. 10[) it follows that S mmo , m ^ to,q and A m 5 mm are bounded and 



have bounded inverse. 



□ 



Remark 5.3. In our previous consideration we assumed implicitly that our se- 
quences {AJ™°} fe z are different from JA™} , i.e. A™ ^ A' fc ™° for all k,k <E 

Z, ?7i, mo <E {1, . . . , , n} in particular A™ 7^ A™ for every k <E Z, m G {1, . . . , n} . 
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Now let us allow A™ = A™ for some indices k € I C Z. Note that in this 
case the operators S mmo , mo y^ rn, are still well-defined and the operator A m S ram 
can be well-defined as well if we define its components as limits of correspondent 
components when A™ — > A™, k G I. This means that for k € I all non-diagonal 
elements of the fc-th line of A m S mm equal and the diagonal elements equal 1. 
Besides, S mmo , mo ^ m and K m S mm remains bounded and with a bounded inverse 
£ 2 — > (-2 operators since formulas (|5.9I) (|5.10l) remain true also when A™ = A™. 
Finally if we consider the dependence A m S mm of the sequence {A™} feeZ , one can 
easy prove that 

A m S mm ({A j. }) — > A m S mm ({X k }) 

as 

^|AT-A™| 2 ^0. 

k 

In other words this means that operators A m S mm and, as a consequence also its 



inverse operators (A m S mm ) , depend continuously of sequence {A™} fceZ on the 



set 

2 



j{An:E| A ™- A £ 



< oo 



Now we are ready to prove our main results on the spectral assignment. 

Theorem 5.4. Let /i 1; [i2, ■ ■ ., /J, n be different nonzero complex numbers and Z\, 
z 2 , •••, z n be nonzero n-dimensional linear independent vectors. Denote 

A™ = In \n m | + i (Arg \x m + 2irk) , m = 1, . . . , n, k S Z. 
Let us consider an arbitrary sequence of different complex numbers {A™} 



fees 
m— l,...,n 

such that 

2 

< oo, m = 1, . . . , n. 



Ek- A ™ 



/,- 



Then there is a small enough e > such for any sequence of nonzero vectors 
{<i™} kel . satisfying 

m— l,...,n 

J2\\d7 -z m \\ 2 <e, m=l,...,n 
k 

one can choose matrices A_i, A2(0) , A^{9) such that for the system (jl.ip , with these 
matrices, the following two conditions hold: 

i) all the numbers {A™} are roots of the characteristic equation det A(A™) = 

0, k € Z, m = 1, . . . , n, 
ii) d™ are right degenerating vectors for A(A™) : rf™*A(A™) = 0, m = 
1, . . . , n; fc G Z. 

Such a choice is unique if we put the following additional condition on matrix Aj,{6) : 



(C) J A 3 (9)d8 = Q. 
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l:S 



Proof. First we denote by A^i the matrix uniquely denned by the relations: 

and denote by {yj} =1 n the bi-orthogonal basis with respect to {zj} =1 in 

C™. Then the corresponding operator A, generated by the matrix A-\, has eigen- 
values A™, m = 1, . . . ,n, k £ Z; A o = and possesses the Riesz basis of eigenvec- 
tors of {^™}fe e z U l^?} — i ■ ^ n ^ ne case wnen an Mm 7^ 1j m = lj.-.jTi, the 
corresponding eigenvectors of .A are 

[ £JT = ( JfeS V m = l,...,n; ft G Z, I^T = A m $T 
(5.11) ^ 



e Afc "j/r, 



In the case /xi = 1, the vector 



' % A -^M, j = l,...,n; I$ = 0. 



^ " 1, e yj 



.7 = 1, 



is a generalized eigenvector of A corresponding to Ao and the other £>™, iffi are 
given by formula (|5.1ip . Let us show that there is a choice of a bounded operator 
To : Xj —> C" (or equivalently a choice of matrices A2(9) , A3 (9)) such that 

(5.12) X 1 ^ ea(A + B Q Po), 
or equivalently det A(A m °) = 0, and 

(5.13) d™°*A(A™°)=0, m = l,...,n; k E Z. 
We represent vectors d™° in the basis {zj} =1 n , namely 



"fen 



E< 

771 — 1 



■mmo*"! 



fco G Z; to, mo = 1) 



With these notations the condition 

E H<C ~ z ™l| 2 < °°) m=l,...,n 
feo 
implies 

E lOmmo I 2 < °°> m ^ m o; E I"™™ - if < OO, m = 1, 

ko feo 

and these sums tend to zero as 

(5.14) J2\K-z m \\ 2 ^0- 



Therefore, under condition (|5.14|) operators D mi m = 1, . . . ,n tend to block diag- 
onal operators 

[S ml ... ... 



D, 



... A m S ; 



m^mm 
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which have a bounded inverse due to Lemma 15.21 This means that for a small 



enough e > 0, the inequality J2 IM™ — z m\\ < £ implies that operators D m , m = 

k 

1, . . . , n have bounded inverse. If numbers {A™} are different from {A™} the later 
fact yields the existence of a bounded operator Vo : Xa — > C™ for which the 
relations (|5.12j) . (|5.13[) are satisfied. Besides, this operator is unique if we require 
additionally: Vo$j = 0, j = 1, . . . ,n, which is equivalent to condition (C). If we 

allow coincidence A™ = A™ for some indices {k, to} £ I one needs to use continuous 
dependence of operators T> on the sequence {A™} (see Remark to Lemma |5.2|) . 
We approximate {A'H by {A'[ n } (E|AT~A™| 2 -► 0) such that {A'H ^ A™. 

k 

Since the conditions (|5. 12|) . (|5.13[) are satisfied for operator "Po({A'™}) they are 
also satisfied for "Po({A™}). This completes the proof □ 

~ 2 
Lemma 5.5. Let {A™} be a given sequence such that ^] [A™ — A™| < oo and 

{®mm } be an arbitrary sequence satisfying 

Y l«mmj < °°) m/mo; ^ |a^ m - 1| < oo. 

kg ko 

Then for any e > and to — 1, . . . ,n there is a sequence {dQ! m } satisfying 

V la fco - S fc ° I 2 < e 

and such that the operator D m has a bounded inverse. 

Proof. First, for given {A™}, let us choose Eq > such that D m will be invcrtible 
for 

Y l a ™V> I' < £ o, m ^ m , Y I a ™« - !| 2 < £ o- 

ko ko 

Then, one can find a great enough N such that 

Y l S ™m O | 2<£ 0, Y l a mm - if < £0- 

|fco|>JV |feo|>JV 

Next we consider the sequences {afy m } for which 

(5.15) al,=«l , as|fc |>iV. 



irint,) 

mnrr fninnnnonfo /t> 

mrriQ i 



Our goal is to choose the remaining components o-m mQl l^ol < N m order to sat- 
isfy the requirement of Lemma. Denote rows of matrix D m by (^™ )*, t™ = 
l,...,n; feo G Z and let g™° be the correspondent components of the vector 
q = D niP , i.e. q^ Q ° = (C °Vp, p G IJ mo = l,... 9 n; Ad € Z. The space 
L = £2 x ... x £2 may be written as L = L 1 © L 2 , where 

n times 

L 1 = {g : q™° = 0, m = l,...,n, |fc | > N}, 
and 

^ 2 = U ■ C = °. N < AT, mo = 1, . . . ,n}. 
Let P be the orthogonal projector on Li- Let us observe that the lines (IT )* for 
\ko\ > N do not depend on chosen components a^ mo , |fco| < iV and that if we put 
a mm = dmmai \ko\ < AT, m, too = 1, . . . , n, then the operator D m : L — > L has 
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a bounded inverse. This means that for all sequences {a£° mo } satisfying (|5.15[) we 
have PD m L = L 2 and the invertibility of D m occurs if and only if 

(5.16) D m L D Li. 

Let L 1 be the subspace 

{peL: (£™<>)*p = 0, |fc | >N,m = l,...,n} 

of dimension (2N + l)n. Denote by £i , j = 1, . . . , n; fc| < iV a basis of L 1 . 
With this notations one can see that (|5.16|) is equivalent to the invertibility of 
(27V + l)n x (27V + l)n matrix 



M 



{(C)*4'. N<iV, |fc|<7V, m = l,...,n, ; = l,...,n} 



i.e. det M ^ 0. The components of M are linear functions of chosen a^ mo , |fcg| < 
TV. Therefore, its determinant is a polynomial of these coefficients. Besides, detM 
is not identical zero because matrix D m is invertible if we take cx%° m = ^rnm a , i-c. 
detM / 0. This implies that M is invertible almost everywhere in £, 2 ( N + 1 ) n . This 
fact completes the proof of Lemma because we can choose a^ mo , \ko\ < N, in such 
a way that 

El fco _^ fc o |< £ 
|fe |<JV 

and the operator D m will be invertible. □ 



Remark 5.6. From the proof of Lemma \5.5\ it is easy to see that actually sequence 
a mm ma V differ from a„ mo only for a finite number of components \ko\ < N. 



Due to Lemma 15.51 the formulation of Theorem 15.41 may be generalized in the 
following way. 

Theorem 5.7. Let the sequence {X™}kez, m=i,...,n an d vectors z m , m = l,...,n 

be chosen according the assumptions of Theorem \5.4\ Then for any sequence of 

vectors 

{d£}kez, m=i,...,n satisfying 

•S—^, ^ 2 

2_^\\d™ - z m \\ < oo, m = l, ...,n, 

k 



and for any e > 0, there is a sequence {d™}/c 



£2 



Y,\\dT-dT\\ 2 <e 

k 

such that, for some choice of matrices A_\, ^4 2 ($)j A 3 (9), satisfying J_ 1 A 3 (8)d9 = 
0, the conditions i), ii) of Theorem \5.4\ are verified. Moreover, {d™} may be chosen 
in such a way, that d™ = d™ for all \k\ > N and for some N G N. 

And then we obtain the following result. 

Theorem 5.8. Let the sequences {A™}fc 6 z, m=i,...,n and {d™}kez, m=l,...,n be from 
Theorem \5.7\ Let, in addition, the complex numbers A°, j = l,...,n be different 
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from each other and different from A™ and let en , j = 1, . . . , n be linear independent 
vectors. Then, for any e > there exist N > 0, a sequence {rf™} feeZ m=1 n '■ 

£ ||C - C|| 2 < e, C = dl\ \k\>N,m = l,...,n 

k 

and a choice of matrices A^i, A2(9), A 3 (9) such that: 

i) all the numbers {A™}fc e z : m=i,...,n U {A^}j=i.... in are roo/Js o/ t/ie character- 
istic equation dct A(A) = 0; 
ii) d™*A(A^) =0, m=l,...,n, fceZ and df A(A°) = 0. 

Proof. Denote by C a (71 x 71) matrix uniquely defined by the equalities: 

dfc = \y;, j = i,..., n , 

and let us put 

ff =df{I-^C)-\ fceZ, ro = l,...,n 

if \l l ^ and 

If =dfC- 1 
for A£ l = 0. 

It is easy to see that the sequences {f™} are also quadratically closed to z m : 

J2\\fl n -z m \\ 2 <oo. 

k 

Therefore, due to Thcorcm l5.7[ for any e > 0, there exist matrices A_i, A 2 (0) , A 3 (9) 
(J. A 3 (8)d9 = 0), a number N > and a sequence of vectors {/™}feez, m=i,...,n : 

E n/* m - /rf < ^, fr = Tr,\k\>N 

k 

such that 

firZ(X?)=0, kGZ,m = l,...,n; 
where M = sup{||I - ^C||, A£ l ^ 0; ||C||} and 

/o (-0 

Xc xe A 2 (9)d9 + c X0 A 3 (9)d9. 

Now, let us put 

A(A) = (/-ic)A(A), A^0. 

One can note that A(0) = and the function Ai(A) = jA(A) may be extended to 
zero by the formula 

f° ~ Z" e Ae - 1 - 

Ai(0) = J-A_!+ / A 2 (fl)d0+ / lim A 3 {9)d9 

7-i J-i*->o A 

Then, one can define 

A(0) = -CA^O). 

Let us observe that A (A) can be written as 

/o /-O 

Ac Ae l 2 ((9)d<9 + / c xe A 3 {6)d9 -C + e _A CA_i 
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r /.0 



c xe CA 2 (8)d8+ c xs CA 3 ( T )dTde, 

/0 /-O 

Ac A <%(6>)d6> + / c X0 A 3 (6)d0 



'-i J-i 

for ^2(0) and ^3(0) given by 

a 2 (0) = 2 2 {e)-{d + i)c-dCA_ x , 

A 3 {B) = M0)+ CA 3 (T)dT - C - CA-i 

It remains to note that, with this choice of matrices, the conditions i), ii) are 
satisfied. Indeed: 

dot A(A) = det(J - -C) det A(A), 
A 

so all the numbers {A™} and {A}° are roots of the characteristic equation and 

df*A(Xf) = rA(A?), k e Z, m = 1, . . . , n, 
= 0, 
dfA(A«) = 0, j = l,...,m 

and finally, for all m = 1, . . . , n, we have 

Y.\\ d T-dT\\ 2 <Y, M2 \\f™-fi n \\< £ - 

k k,m 

This completes the proof. □ 



6. Conclusion 

We give here some conditions on sets of complex numbers {A} and n- vectors {d} 
such that they form a spectral set for a neutral type systems. This is a first etap 
for solving vector moment problems using the exact controllability properties of a 
neutral type system related to the given moment problem. 
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